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STRAIN-HARDENING MATERIALS 
by Alexander Mendelson 
Lewis Research Center 

SUMMARY 

A simple and straightforward procedure is presented for solving the elastic -plastic 
problem for the torsion of a solid prismatic bar made of a strain -hardening material. 
The procedure is based on the method of ’’successive elastic solutions” or successive 
approximations. For circular cross sections, the problem reduces to the solution of a 
nonlinear algebraic equation. For the case of a circular cross section with linear strain 
hardening, the solution is obtained in closed form. Results are presented for bars of 
rectangular and circular cross sections with linear strain hardening. 


INTRODUCTION 

Although the problem of elastic torsion of prismatic bars is one of the classical 
problems of mechanics and has received extensive treatment, the corresponding problem 
of elastic -plastic torsion has as yet not been satisfactorily attacked. 

The usual solutions to the plastic torsion problem assume perfectly plastic materials 
and are generally limited to cases in which the complete cross section is plastic, because 
the elastic -plastic boundary is considered difficult to find. Furthermore, for completely 
plastic sections of bars made of perfectly plastic materials, the ingenious sand-hill 
analogy of Nadai (ref. 1) enables one to obtain solutions experimentally. Examples of 
solutions of this type can be found in references 1 to 5. An example of a solution for a 
strain-hardening material can be found in reference 6 for a special type of stress -strain 
law. 

This report presents a relatively simple solution to the elastic -plastic torsion 
problem for strain-hardening materials. Use is made of the method of successive elastic 
solutions for successive approximations which have been applied so successfully to many 
other types of problems as described in references 7 to 10, among others. For the 



practical and important case of a bar of circular cross section, the problem reduces to 
the solution of a nonlinear algebraic equation. For linear strain hardening, the solution 
is obtained in closed form . 


METHOD OF ANALYSIS 
Basic Equations 

The analysis begins with the Saint-Venant semi-inverse method as for the elastic 
case. Consider a prismatic bar subjected to a twisting couple as shown in figure 1. One 



end of the bar is assumed fixed against rotation but not against warping, and, at the other 
end, a couple M with moment along the z-axis is applied. Following Saint-Venant, it 
is assumed that the displacements in the x, y, and z directions, respectively, are 
given by 


u = - yz a 

V = XZQ! S (1) 


w = w (x, y, a) j 


where a is the angle of twist per unit length. (All symbols are defined in appendix A.) 
The warping function w is directly proportional to the angle of twist per unit length a 
in the elastic case. For the plasticity problem this is, in general, no longer true. 
Substituting equations (1) into the usual strain -displacement relations results in 
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( 2 ) 


yxz = - «y + 


3w 

9x 


y = ax + 

r yz 


3w 

3y 


all the other strains being zero. If equations (2) are now substituted into Saint- Venant's 
six compatibility equations, four of them are identically satisfied and the remaining two 
become 


which gives 


3 y 


yz 


3y 


xz 


0X \ 

3x 

3y J 


3y 

y z 

3 y 

! A XZ 

3y l 

3x 

9 y , 


= 0 



= constant 
3y 


( 3 ) 


(4) 


From equations (2), it follows that the constant appearing in equation (4) must equal 2a. 
Hence, 



9x 



2 a 


Equation (5) is the compatibility equation for this problem. 
The only equilibrium equation not satisfied identically is 


T 3t 

xz + y z = o 

0x 3y 


The stress -strain relations can be written 


( 5 ) 


( 6 ) 


3 


1 p 

y = — T + y c 
r xz q XZ XZ 


> 


y = 1 r p + yP 

y z g y z yzj 


( 7 ) 


where y£ and y£ are the total accumulated plastic shear strains. Substituting equa- 
xz yz 

tion (7) into equation (5) yields the compatibility equation in terms of the stresses 


where 


9t „ 9t 

— — = 2Ga + g(x, y) 

9x 9y 


g(x, y) = G 


'd y 


P 

xz 


3y 



( 8 ) 


0) 


The stress function is now introduced such that 

T = d Jl T = - d Jl 
xz 9y ^ 9x 

Then the equilibrium equation (6) is identically satisfied, and the compatibility equa- 
tion (8) becomes 


(10) 


V 2 (^ = §?£= - 2G« - g(x,y) (n) 

_ 2 _ 2 
9x 9y 

For the elastic problem g(x, y) is equal to zero. 

The boundary condition for this problem can be obtained directly from the condition 
that the lateral surface is force free. As shown in any standard text this becomes simply 

cp = 0 ( 12 ) 

on the boundary, provided that the section is simply connected. For multiply -connected 
cross sections, the stress function <p is equal to a different constant on each boundary. 
We will concern ourselves here only with simply connected sections, that is, solid bars. 
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From equation (10) it follows that the resultant shear stress is equal to the gradient 
of <p ; that is, 



The resultant moment can readily be found by integration to be 

M = 2 II < p dx dy (14) 

The elastic -plastic torsion problem has thus been reduced to solving equation (11), 
subject to boundary condition (12). For the elasticity problem, the right side of 
equation (11) is a constant, and solutions, at least for simple shapes, can readily be 
obtained. For the plasticity problem, the right side of equation (11) is a function of the 
plastic strains and is therefore unknown until the solution is obtained. To solve the 
problem, use is made of the method of successive elastic solutions or successive approx- 
imations. 

Before proceeding with the details of the solution, the following should be noted: 

As shown in reference 11, both the total and incremental theories of plasticity furnish 
the same solution to the torsion problem provided either (1) the cross section is circular, 
or (2) the material is perfectly plastic. It is reasonable to assume, therefore, that this 
will be approximately true for most practical problems. Indeed, it has been shown in 
reference 6 that for the case of a square cross section with strain hardening there is 
little difference between total and incremental theories. In what follows, therefore, use 
will be made of the total or deformation theories of plasticity and the load will be as- 
sumed to be applied in one step. The use of incremental theories does not appreciably 
complicate the problem, and the necessary formulation is given in appendix B for those 
who desire to use it. 

The deformation theory of plasticity using the von Mises yield criterion results in 
the following plastic strain - total strain relations as shown in reference 7. 


(15) 


where 
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( 16 ) 



and the equivalent plastic strain e p is related to the equivalent total strain e t through 
the stress -strain curve of the material and can be written, in general, as 


e 


P 


= f(e t ) 


(17) 


These plasticity relations will now be used for obtaining a solution by successive 
approximations . 


Method of Solution 


It is convenient to introduce the following dimensionless quantities: 


U = M* = 


M 


2Ge Q a 


_ y xz 
r x = — 
x 2e„ 


y 

y *0 


e p s r 

e 0 


T„ = 


xz 
2Gen 


yP 

O E ’"xz 


2e r 


2Ge 0 a 

T 2e 
ze 0 


e t B 


T = 

y 2Ge 


0 


5 ^ v = y - 

a a 




(18) 


where e Q and cr 0 are the yield strain and yield stress, respectively, in the uniaxial 
tensile test, related to each other by ctq = E€q, and a is a characteristic dimension of 
the cross section. 

The system of equations to be solved for a simply connected cross section can now 
be written as follows: 
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(19a) 


g(k,v) 



877 34 


V^U= - /3 - gU.rf) 

U = 0 on boundary 


T x = 


3U 

3y 


r„ = - 


3U 

34 


r x = T x + lx y y = T y + ly 


e t = ^\ yz + r 


2 , 2 
^f' x r y 


e p = F < e t> 


^ = J2 y yP = _B ^ 

X 'x 'y 7 t 

e t e t 


(19b) 

(19c) 

(19d) 

(19e) 

(19f) 

(19g) 


The successive approximation method proceeds in the following manner. The 
plastic strains are assumed to be zero everywhere. Equation (19b) is solved by any 
available method. The stresses, the total strains, and equivalent total and plastic 
strains are computed by means of equations (19c) to (19f) with the help of the stress- 
strain curve. If, at any point in the cross section, the equivalent plastic strain as 
computed from equation (19f) is negative, this point is in the elastic region, and the 
plastic strains at this point are set equal to zero. Otherwise, new approximations to the 
plastic strains are calculated by means of equations (19g). One then returns to equa- 
tion (19a), and the process is repeated until the difference between two successive 
iterations differs by less than some preassigned value. For all the cases treated here- 
in, convergence was obtained very rapidly. A discussion of the convergence of this 
technique for other types of problems can be found in references 7 and 8. The method 
will be illustrated for bars of rectangular and circular cross sections. For a circular 
cross section with linear strain hardening, the solution can be obtained in closed form. 
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Bar With Rectangular Cross Section 


Equation (19b) is the well-known Poisson equation of mathematical physics which can 
be solved by a variety of methods. Herein the finite -difference method is used because it 
is simple and straightforward and is readily programmed for a digital computer. 

Consider a bar of rectangular cross section as shown in figure 2. Note first that, 
because of symmetry, only one quadrant of the section need be considered. For a 
square cross section, the diagonals are also lines of symmetry, and only one octant is 
used. The quadrant is divided into a grid of nxm points as shown in figure 3. At each 
of the grid points such as the point designated by (i, j) equation (19b) is written in finite 
difference form 


U, 


i+lj 


U. , - + U. . i + U. - n 

l-l, j i,j-l i,l+l 


4U. 

i, 


-03+ §ij) h ' 


(20) 


where h is the constant grid spacing divided by the half width of the plate a and is the 
same in both the x and y directions, and 


T,l 


_1_ 

2h 


yP ... -yP 
x, i,]+l ; 


x,i, j-1 " y y,i+l, j + 


yP . 1 . 

y y, i-i,] 


(21) 


An equation such as equation (20) can be written for each of the mn grid points 
resulting in a set of mn simultaneous linear equations for the unknown values of U 
at each of the points. Actually, the number of equations to be solved is (m - 1) (n - 1) 
rather than mn, since the boundary conditions require U to be zero at the upper and 




b 


a 


-St= x 


Figure 2. - Rectangular cross section. 



Figure 3. - Finite-difference grid for one quadrant of 
rectangle. 
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right hand boundaries of the quadrant. For a square section, only the points on the diag- 
onal and to the right of the diagonal (or left) need be considered, resulting in l/2n(n - 1) 
equations. Along the lower boundary, because of symmetry, equation (20) becomes 


U i + l,l + U i-l, 


+ 2U 


i, 2 


4U i, ! = - 0 3 + gi i)h' 


L 


Along the left boundary, 


U 


i,j-i + u i,j + i +2u 2,j- 4u i,r- (,! + B i.j )h 


at the center, 

4U 2, 1 - 4U 11 = - SlPh 2 


and along the diagonal of a square, 


2 u i + l,i + 2 u i,i-l- 4 V = - fti+g M> h ‘ 


Once the values of U are determined at all the grid points, corresponding values of 
plastic strains are computed by means of equations (19c) to (19g). The g^ j are then 
recomputed, and equations (19b) are solved again. The process is repeated until con- 
vergence is obtained. An example using this technique is described in the section 
RESULTS AND DISCUSSION. 


Bar With Circular Cross Section 

For a bar with a circular cross section, the solution is greatly simplified. In 
particular, for the case of linear strain hardening, a closed-form solution can be ob- 
tained. In polar coordinates, the displacements are 
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( 22 ) 



and the only nonzero strain is 


y 


0 z 


= ar 


The stress -strain relation can therefore be written 


1 

G 


6 z 


= ar - y 


P 

6z 


(23) 


(24) 


The von Mises equivalent stress reduces to 

CT e = V3 = y[S Gar - yfiGy^ (25) 

and the equivalent plastic strain is 


Hence, 



P 

6z 


(26) 


Let 


o e = yJdGar - 3Ge p 


(27) 



— = (3 e 
e 0 



(28) 


where a is the radius of the bar. Then equation (27) can be written in dimensionless 
form as 


& 

2 


PP ~ - e p 
2 p 


(29) 
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The stress -strain curve equation can be written in dimensionless form as 

s e = ffep) (30) 

and combining with equation (29) results in 

2f(e p ) = VTfr) - 3e p (31) 

which can be solved iteratively for e p . 

Equation (31) is valid only in the plastic region. Let this region extend between 

p = P c and p = 1. To determine the position of the elastic -plastic boundary p c , let 

e =0 when a = a n or when 
p e 0 


s e = ! + P 


Hence, from equation (29), 


_ 2(1 + p) 

C 


which depends only on Poisson’s ratio p and the yield strain but not on the stress- 
strain curve. The value of /3 at which plastic flow just starts is found from equa- 
tion (32) by setting p £ equal to 1. Thus, the critical value of /3 will be 


B =20±P) 

C y/2 

and the critical angle of twist per unit length will be 




„ - 2 (1 + M ) £ o 

c a c " 


e r 
rtZ 

a 


or 


n, r 1 CT ° 
O' Or = — — 

c v/3 a 


(33) 


(34) 


In summary, the strain-hardening solution is found as follows: The elastic -plastic 
boundary p c is first determined from equation (32). The stress and strain in the 
elastic region for p < p c are then computed from equations (23) and (24) with y^ z set 
equal to zero. In the plastic region p > p c equation (31) is usually solved by an 
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iterative method. Then, yjj can be computed from equation (26), and the shear stress 
from equation (24) . Once the shear stress is known throughout the section, the torque 
can be computed by integration. 

Let us now consider the case of linear strain-hardening. Equation (30) for the 
stress -strain curve, can be written as follows: 

S e = (1 + M)(l + e p) (35) 


(where the strain-hardening parameter m is the ratio of the slope of the linear 
hardening curve to the slope of the elastic curve). Hence, equation (31) becomes 


or 


2(1 + M ) 



m 

1 - m 



V3~fl p - 3e p 


e - 2 (1 + H) 

P , , 2(1 + M)m P_P ' 
1 - m 


(36) 


Note that the critical value of p is obtained when the numerator of equation (36) 
vanishes, which results again in equation (32). 

Once the equivalent plastic strain is known from equation (36), the plastic shear 
strain and the stress are computed from equations (26) and (24)„ Thus, a complete 
solution in closed form is obtained. To compute the torque, define 

T e s W 2Ge o 


(37) 


t 9 =±Pp p^p c 
T e = l P p - e e p - p c 


Then, 


M* = M/2Ge 0 a 3 = 2i t J r 0 p 2 dp 


and substituting 


12 


J 



results in 


where 


M*-f- V^[j A(l - 4). i B(l - p3)] 

V3“/3 


A = 


3 + 


2(1 + p)m 
1 - m 




( 38 ) 


B = - 2 (! + M) J 

3 + 2(1 + 

1 - m 

Note that for P c = 1 (no plastic flow) the torque reduces to the elastic torque as given 
in standard texts. For a perfectly plastic material 

m = 0 A = -£- B = - — (1 + p) 

3 

and 


or 


M* = 1 -ip 3 ) 

3^3 l 4 c / 


M = 


27ra 3 OQ 

3^3 


1 


12 yfi e 



3-| 


and the classical solution as given in reference 2 is recovered. 


(39) 


RESULTS AND DISCUSSION 

Calculations were performed by the above techniques for two illustrative cases: The 
case of a bar of square cross section and the case of a bar of circular cross section. A 
value of 0. 3 was used for Poisson's ratio for all calculations. For the square cross 
section, the finite -difference solution previously described was used with an 11 by 11 
grid as shown in figure 3. As a check on the accuracy of the finite -difference formula- 
tion, the results of the elastic solution which were obtained from the first iteration are 
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shown in the following table. These results are compared with those of reference 12. 



Reference 12 

Finite 

differences 

r „ /Gara 
max 7 

1. 351 

1. 343 

M/Gaa^ 

2.250 

2.244 

<Pl j/Gaa 2 

. 589 

.589 


It is seen that the solution with this many grid points is sufficiently accurate. 

The elastic -plastic calculations were performed assuming linear strain hardening. 
For linear strain hardening, equation (19f) can be written as 


e 


P 


e t - | (! + M) 
o 

1 + 1 (1 + 

3 1 - m 


(40) 


where the strain-hardening parameter m is the ratio of the slope of the linear strain- 
hardening curve to the slope of the elastic curve. For the perfectly plastic case, m 
is equal to 0 and, for the elastic case, m is equal to 1. Figures 4 to 7 show the 
effects of the strain -hardening parameter and the angle of twist on the maximum stress, 



Figure 4. - Variation of dimensionless maximum shear stress with 
dimensionless angle of twist per unit length for several values of 
strain-hardening parameter for square cross section. 
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Dimensionless angle of twist per unit length, p 

Figure 5. - Variation of dimensionless maximum plastic 
shear strain with dimensionless angle of twist per 
unit length for several values of strain-hardening 
parameter for square cross section. 



Dimensionless angle of twist per unit length, p 


Figure 6. - Variation of dimensionless moment with dimensionless 
angle of twist per unit length for several values of the strain- 
hardening parameter for square cross section. 
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I I II 


TABLE I. - SUMMARY OF RESULTS FOR TORSION OF 
SQUARE PRISMATIC BAR 


Linear 

Dimensionless 

Dimensionless 

Maximum 

Maximum 

strain- 

angle of twist 

torque, 

dimensionless 

dimensionless 

hardening 

per unit 

M* 

shear stress, 

strain, 

parameter, 

length, 


T max 

■)P 

'max 

m 

P 





0 

2 

1.786 

0. 751 

0. 820 


3 

1.918 



1. 824 


4 

1.955 



2. 851 


K 

1.977 



3.959 

0.05 

2 

1.813 

0. 785 

0. 758 


3 

1.997 

. 825 

1.623 


4 

2.094 

. 862 

2.434 


5 

2. 166 

. 899 

3.240 


6 

2.228 

. 934 

4.003 

0. 10 

2 

1. 838 

0. 818 

0.7015 


3 

2.073 

.893 

1.478 


4 

2.223 

. 963 

2.209 


5 

2.347 

1.032 

2.919 


6 

2.465 

1.099 

3.618 

0.20 

2 

1.890 

0.881 

0. 600 


3 

2.220 

1.022 

1.250 


4 

2.471 

1. 156 

1. 870 


5 

2.717 

1.290 

2.488 


6 

2.966 

1.426 

3. 116 
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Figure 7. - Plastic zone boundaries in quadrant of square 
cross section as function of dimensionless angle of twist 
per unit length for strain-hardening parameter, 0, 1, 



the maximum plastic strain, the torque, and the size of the plastic zone. The results 
are also summarized in table I. 

The results of the computations for a circular cross section are shown in figures 8 
to 10. Linear strain hardening was assumed, and equations (32), (23), (24), (36), (26), 
and (38) were used. Figure 8 shows the elastic-plastic boundary as a function of (3. 
Figure 9 shows the effect of the strain-hardening parameter on the shear stress for 
j3 = 5. 0, and figure 10 shows the relation between these parameters and the torque. The 
corresponding quantities for the other values of m and j3 can easily be computed. 

It is seen from the preceding discussion that solutions to the elastic -plastic torsion 
problem can be obtained in a relatively simple and straightforward manner. It is not 
necessary to assume that there is perfect plasticity or that the complete section has 



Dimensionless elastic-plastic boundary, p c 

Figure 8. - Variation of dimensionless elastic-plastic boundary 
with dimensionless angle of twist per unit length for circular 
cross section. 
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Figure 9. - Variation of dimensionless shear stress with 
dimensionless radius in circular cross section for 
various values of strain-hardening parameter. 
Dimensionless angle of twist per unit length, 5.0. 



Figure 10. - Variation of dimensionless moment with 
dimensionless angle of twist per unit length for 
various values of strain-hardening parameter 
for circular cross section. 

yielded. The elastic -plastic boundary is found automatically as part of the calculation 
and does not pose any special difficulties in its determination. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, October 4, 1967, 
129-03-08-04-22. 



APPENDIX A 


SYMBOLS 


A, B 

constants 

e 

strain 

a 

characteristic dimension 

M 

Poisson’s ratio 


of cross section 


dimensionless coor- 

e 

"dimensionless" strain 


dinates, eq. (18) 

G 

shear modulus 

P 

dimensionless radius, 

g 

plastic strain function 


eq. (28) 

M 

torque acting on cross 

a 

stress 


section 

T 

resultant shear stress 

M* 

dimensionless torque, 

T xz’ T yz’ t 8z 

shear stresses 


eq. (18) 

T x’ T y’ T e 

dimensionless shear 

m 

linear strain-hardening 


stresses, eq. (18) 


parameter 

<P 

stress function 

r 

S p 

radial coordinate of 
circular cross section 

dimensionless equi- 

Subscripts: 

c 

critical 

e 

valent stress 

e 

equivalent 

U 

dimensionless stress 

i,j 

station index 


function, eq. (18) 

P 

equivalent plastic 

U, V, w 

displacements 

r 

radial 

x,y,z 

coordinates 

t 

equivalent total 

a 

angle of twist per unit 

e 

tangential 


length 

0 

yield 

P 

dimensionless angle of 
twist per unit length, 

Superscripts: 



eq. (18) 

P 

plastic 


r xz’ r yz ,y 0z en S ineerin g shear strains 

y ,y ,y "dimensionless" shear 

r x ,r y ,r z 

strains, eq. (18) 



APPENDIX B 


INCREMENTAL PLASTICITY EQUATIONS 


If the incremental or flow theories of plasticity are to be used, the angle of twist 
per unit length a must be increased in small steps and the same calculation as pre- 
viously described performed after each step. The stress -strain relations (eq. (7)) are 
written as follows: 


y = — t + y p + Ay P 
'xz G xz 'xz 


xz 


y = — t + y P + Ay P 
r yz Q yz 'yz 'yz 


(41) 


J 


where y p and y p are the accumulated plastic shear strains up to, but not including, 
xz y z « Y) 

the present increment of load, and A y£ z and Ay^ z are the increments of plastic shear 
strains due to the current increment of load. The equation for the stress function cp 
now becomes 


V c p = -2G(ai + Aa) - g(x,y) - Ag 


(42) 


where g(x, y) is as previously defined, and 


Ag = G 


0/Ay 


xz 


0 Ay 


n, 


9y 


0X 


(43) 


g(x,y) is known, and Ag is to be determined. Equations (15) to (17) are replaced by 
the following: 


Ay; 


P - A£ P .,P 


xz 


-t 


Ae 


yz 


r. 


yz 


(44) 


Ay p = — E y p 


yz 


j 
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where 



( 45 ) 


(46) 


(47) 


(48) 


where a g , appearing in equation (47), is the value at the end of the previous increment 
of load and is known. Similarly, da /de is the slope obtained from the uniaxial stress- 

P 

strain curve at the end of the previous increment of load and is known. Equation (47) 
is approximate except for linear strain hardening when it is exact and can be written as 


Ae 


P 


2 1 + ji 

3 E 


a„ 


l + - (i + a) 
3 


m 

1 - m 


(49) 


The problem can now be solved for each increment of load by successive approxima- 
tion as described in detail in references 8 and 9. For a given increment Aa, values 

of AyP and A y? (such as zero) are assumed. Equation (42) is now solved for cp 

xz yz ' 

since Ag (as well as g(x, y)) is known. The stresses can be computed from equa- 
tions (10), the strains from equation (41), the modified total strains from equations (46), 

e. from equation (45), Ae from equation (47), and new values of Ay* 5 and Ay* 5 
* p xz yz 

from equation (44). Then, Ag is recomputed and the process repeated until convergence 
is obtained. 

By the above process, the plastic shear increments can readily be obtained for each 
increment of load. However, as previously pointed out, such an incremental process is 
not necessary for circular cross sections or any cross sections of perfectly plastic 
materials. 
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